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Abstract. We study the question weather weak solutions to a class of active scalar equations, with 
the drift velocity and the active scalar related via a Fourier multiplier of order zero, are unique. Due 
to some recent results we cannot expect weak solutions to be unique without additional conditions. 
We analyze the case of some integrability conditions on the gradient of the solutions. The condition 
is weaker than simply imposing £ L°° . Lastly, we consider the inviscid limit for the studied 
class of equations. 



1. Introduction 

Active scalar evolution equations have been a topic of considerable study in recent years, in part 
because they arise in many physical models. In particular, such equations are prevalent in fluid 
dynamics. This paper is concerned with a class of diffusive transport equations for an unknown 
scalar field 6{x,t), of the form 

dte + {u ■ v)e + K{-Ay9 = o, (dde) 

where k > 0, and their non-diffusive version 

dt9 + (n • V)e = 0, (DE) 

For both equations we assume t > and x € M". The divergence free velocity field u{x,t) 

divu = (1) 

is obtained from 6 via C alder on- Zygmund singular integral operators, defined in the frequency 
space by a multiplier, whose symbol is of order zero (see ^H] for a precise definition) 

m = mm. (2) 

Examples covered by the above class of equations arise in a variety of physical models. To be more 
specific we present a few examples of such. 

Example 1.1. The 2D surface quasi-geostrophic equation models the evolution of buoyancy (or 
potential temperature) 6 on the 2D horizontal boundaries of general 3D quasi-geostrophic equations 
(see [1] for more details). It is given by 

dte + {u-v)e + {-A)"9 = o, 

with divM = 0. The drift velocity u is recovered from the buoyancy 9 through the multiplier 
operator 
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Example 1.2. The 2D porous media equation governs the motion of the density p of an incompress- 
ible fluid: 

dtp + (n • V)p = 0, 

with div n = 0, where u and p are related via Darcy's law: 

u = -Vp-{0,p). 

For simplicity we have set all physical constants equal to 1. Eliminating the pressure from the 
equation, one obtains u = T[p], where r[-] is a Fourier multiplier operator with symbol 

icr'(eie2,-(6)'). 

In the present paper, we shall be concerned with weak solutions. We say that a pair {6,u) G 
Lf^ciW^ X R) is a weak solution of the equation (|DDE]l on M if for every ^ € C^{W x R) the 
following holds 

[ [ e{dt(j) + (n • V)^ + K{-Ay^) dxdt = 0. 
Jr jr" 

The primary motivation for studying weak solutions in the context of hydrodynamics comes from 
the desire to understand the laws of turbulence in the limit of infinite Reynolds number. For the 
above mentioned examples of active scalar equations, existence of weak solutions has been shown 
for instance in |17] and [3], respectively. For the general dissipative case addressed in this paper, 
we give a proof of existence of weak solutions in the appendix. Proofs of existence for the non- 
dissipative equations take into account the special structure of the non-linearity (the actual symbols 
of the multipliers and related to them commutator estimates) thus we refer the interested reader 
to papers devoted to the analysis of a specific problem. 

The classical way to prove uniqueness of such solutions relays on applying energy methods to an 
equation governing the difference of two solutions, 9 = 9i — 92- The most direct condition which 
guarantees the uniqueness of 9 is to require, for instance, that V6'i € L^{0,T, L°°{W^)). We then 
obtain the following immediate estimate 

< C'||V6'i||^i((o,T),L°°(IR"))Pllioo((o,T),L2(Mn)), (3) 

which is the core of the energy approach. Given this estimate, uniqueness of solutions follows from 
an application of Gronwall's inequality. Note that we impose high regularity on the solution (and it 
is not clear that it follows from energy methods). Moreover, for equations in 2D, the classical result 
gives local existence of solutions 9 with data in Sobolev space with s > 2. Even for solutions in 
H^, we have V9 G H^. This, however, is not enough to conclude that S/9 G L°° and motivates the 
choice of BMO (see next section for the definition) rather than L°°. 

Results of Cordoba, Faraco, Gancedo [9] and Shvydkoy [18] imply that it is not enough to require 
only that {9,u) € L°°(R" x R) (obviously contained in Lf^^iR."- x R)) in order to obtain unique 
weak solutions. Therefore, one needs additional conditions in order to ensure uniqueness of such 
solutions. 

In our paper we consider a situation where \/9 = g + b with g € L^((0, T), L°°(R")) and b G 
LP({0,T), BMO) such that supp b has finite measure for a.e. t G (0,T). More precisely, in Section 
3, we prove that under such condition, where p is related to the strength of dissipation 7, weak 
solutions to (IDDEp are unique. We prove an analogous result in the case of the inviscid problem 
(IDEp . however here we need to impose additional integrability of solutions (since it is no longer 
provided by dissipation). In Section 4, we consider the inviscid limit of (IDDEp . Our proofs follow 
techniques introduced in [l6] . 



9-\/9i9 dxdt 
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2. Preliminaries 

We begin by recalling the definition of the space BMO. Given a function / G L\^^{W^) and a 
cube Q, let fq denote the average of / on Q 

Define the sharp maximal function by 

M#/(x) = sup^ / \f{x)-fQ\dx, 

where the supremum is taken over all cubes Q containing x. We say that / has bounded mean 
oscillation if the function M^^f is bounded, that is 

BMO = {/ G Ll^O^"") : M*f G L°^(R")}. 

We define a norm on BAI O by 

\\f\\BMO = \\M*f\\L^. 

Clearly, L°° C BMO. Recall however that C BMO due to the fact that BMO contains also, 
for instance, unbounded functions with logarithmic singularities, like 

— In \x\ for < Ixl < 1 



for \x\ > 1. 



As shown by our results, the presence of such logarithmic (or in general BMO) singularities, 
localized to a set of finite measure, does not affect the uniqueness of weak solutions. 
In the proofs we relay on the following remark. 

Remark 2.1. Observe that / G BMO is by definition equivalent to M^^f G L°°. Because of the 
above comment, the reverse implication is not true. However, Fefferman and Stein in [12] proved 
that M*/ G implies / G L^, if 1 < p < oo. Moreover, the John-Nirenberg inequality yields that 
if / is in BMO, then / is locally in for any 1 < p < oo. 

It is not obvious that one can consider a dual pairing of L^-functions with functions in BMO 
because (BMO)* = IH} C (the Hardy space). This problem was considered by P.B. Mucha and 
the author in [16]. We recall the result since it will be the main tool used in the proofs. 

Theorem 2.2. Let f G BMO, such that suppf has finite measure and g G L^(M") n L°°(M"). 
Then the following inequality is true 



fg dx 



< CWfllBMoWgh^ [1 + I In llgll^i I + ln(l + IbllL-)] • (4) 



3. Uniqueness of weak solutions 

We first state and prove the result for diffusive equations. 

Theorem 3.1. Let 0{x,t) be a weak solution of the equation WDB^ such that 

6 G L°°((0,T),l2(R")) nL2((0,r),i/^(M")). 

Assume in addition VO = g + b where g G Li((0, T), L°°(M")) and b G L''{{0,T), BMO) for some 
s> where q > 2 (s = oo if q = 2) is such that 6 G L'?((0, T), LP(M")) for some p G (2, ;^), 
with supp b of finite measure for a.e. t G (0,T). Then {6,u) is unique. 
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Proof. Let 9i and 62 be two solutions of (jPDEp and mi,M2 corresponding drift velocities. The 
difference 6 = 61 — 62 satisfies the following equation 

Ot + {ui - U2) ■ vOi + u2-ve + i^i-Aye = o. (5) 

Let us denote A = (— A)^/^. Multiplying by 9 and integrating in space, using the divergence-free 
condition on U2, we obtain 

~[ \e\^dx + vf \K"'e\^dx+f {U1-U2) -s/eiO dx = 0. (6) 

2 dt J^n J^n J^n 

Set Q = (ui — U2)9 and /3 = V6'i. Notice that \\6\\l2 and \\u\\i^2 are both C([0,T]) upon ([6]) and 
assumptions on 6i and 02 • We split a = Om + ctr, where \am\ = min(|a|,m) for some m > 1. 
Let f3 = /3i + /32, where f3i corresponds to the part g and /32 corresponds to the part b of the 
gradient. Theorem 12.21 yields 



\a(3\ dx 



< ll/3i|lL°°||a;||ii + C\\(32\\BMo\\am\\Li (1 + |ln||am||Li| + ln(l + m)) + / 



(7) 



arf32\ dx. 



Notice that ||am||Li ^ IIq^IIli — II^IIl^ H^H^a < Ci|10||^2- Because of continuity of \\0\\l2 and 11^*11^2 
in time, we can pick a time Tq > such that for t £ [0, Tq] 

||am||ii|ln||am||iil < Ci 116*1112 1 In Ci 116*11^2 1, (8) 

thus we get 

\af3\ dx 



< C{\\Mloo + \\I32\\bmo)Ci\\9\\1, (1 + llnCill^llial +ln(l + m)) + [ 



(9) 



|Qr/32| dx. 



Inequality ([9]) combined with ([6]), after dropping the positive dissipation term, give 

1 d ^ i\a\\2 

2crd^^^"^"^^ 

< C{\\Ml^ + \mBMo)Ci\\e\\l2 (1 + [InCill^llial +ln(l + m)) + [ |a,/32| dx. 



(10) 



Denote 



/(t) = 2CiC(||/3i||L^ + ||/32||BA/o), r{t) = Cif |a,/32| dx, x{t) = Ci\\e\\l,. 



Notice that the assumptions on V6 guarantee that f{t) is integrable on bounded time intervals. 
We obtain the following differential inequality 

x'{t) < f{t)x{t){l + I lnx(t)| + ln(l + m)) + r{t), 

x{0) = 0. ^ ^ 

In order to find a good estimate on x{t) we introduce the following equation 

y'{t) = f{t)y{t){l + I lny{t)\ + ln(l + m)) + r(t), ^^^^ 
y(0) = 

for some m large enough. From the Osgood existence theorem we know that there exists a unique 
local solution to ()12p . In addition, the right hand-side of (jl2p guarantees that y{t) is increasing. 
This implies that the solution of ()12p dominates x{t), that is 

< x{t) < y{t) < 1 fortG[0,ro]. 
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Therefore, we investigate the behavior of solutions to (jl2|) . By Gronwall's inequahty we obtain 
y{t) < ^ exp (^j^ /(s)(l + I hiy(s)| + ln(l + m)) ds^ 

+ j r{s)exp(^j^ /(r)(|lny(T)| + 1 + ln(l + m)) dr^ ds 

< ^exp^(l + ln(l + m))^ /(s) ds^ exp f {s)\ In y{s)\ ds^ 

+ exp^^ f {t)\ In y{T) \ dr]^ exp ^(1 + ln(l + m)) ^ f{T)dT^J^ r(s) ds. 

Since 1 > y{t) > 1/m imphes | lny(t)| < Inm we can estimate the right hand-side of ()13p (up to a 
constant) by 

— (1 + m)^o /(s) ds^J^fis) ds ^ ^/o'/W + ^)/o /W / ^(5) ds 

t (14) 
< (2m2)/o '^'(—+ [ r{s) ds) = (2m2)/o /(^) ds f ]_ ^ ^jf f j^^^^j A _ 
Km Jo J \m Jq J^n ) 

In order to complete the proof it suffices to control the part involving . A suitable decay of this 
term is due to the fact that the support of a^. shrinks with growing m. To take advantage of this 
fact additional integrability of a is needed. For the dissipative case this is a consequence of the 
energy estimates and interpolation. More precisely, we have 

aGL'?((0,r),Lf(M")) (15) 

for some 1 < g < oo and 1 < p < ^^^^ given by the usual linear interpolation relation. For e > 
such that p — e > 1, by Holder's inequality, we get 

\arh\ dx < Warhp-AlMiip-,)' i^suppar) - C'll^r ||/32 HbA/O (16) 

Inequality (I16p is true for a.e. t since the support of Or has finite measure and, by Remark 12. 1^ 
functions with bounded mean oscillation are p-integrable on sets of finite measure for all p < oo. 
From Chebyschev's inequality we obtain 



\suppar\ < < . (17) 

ml \ m I 



By Holder's inequality 



1 

l^rllLP-'! = ( / WrT~^ ] < \supp ar\P^P-^ WarWlP (18) 

' suppar 



Combining (jl7l) and (llSp we obtain 



Hence 



\ar(32\dx<m p-^\\a\\j;^p'' 'WPWbmo = m p-' \\a\\lp' \\P\\bmo- (20) 



Therefore 

/ \arP2\ dxds <m f-HI«llL9((0,T),LP(R"))ll/?2|li(7(p-e)/p)'((0T) SMO)- (^1) 

Jr" \\ ' " 
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Hence, if (32 G L'^{{0,T), BMO) for some s such that q' < s < oo, there exists an e > such that 
the above conditions are satisfied and s' = q • Denoting 

M = 11/32 ||l''({o,t),ba/o)) 

estimates p^ - (f2T]) give 

y{t) < C(2m2)/o/W'^- (^^ +m-^^CiM||a||i.((o,T),Lp(R"))) • (22) 

We choose < ti < Tq small enough so that 2 /q' f{s)ds-l<--^^. Then 

y{t) < Km'^^ for < t < ti, (23) 

where is a constant independent of m. Letting m — > oo we get y{t) = for < t < ti which 
implies x{t) = for < t < ti. Therefore 9i = O2 for < t < ti. We can continue this process 
starting at t = ti and obtain uniqueness of 9 for all t € [0,T]. Uniqueness of u follows from 
uniqueness of 0. □ 

In the inviscid case k = 0, the equation does not provide any additional integrability (there is 
no dissipation). Therefore, in order to prove a similar theorem we have to impose this condition 
on our solutions. The following theorem is a corollary from Theorem 13.11 

Theorem 3.2. Let 9{x,t) be a weak solution of the equation liDH^) such that 

6 e L~((0,r),L2(M")) nL'?((0,T),LP(M")) with p,qe (2,00). 

Assume in addition VO = g + b where g € Li((0, T), L~(M")) and b € L'{{Q,T),BMO) for some 
s > (s = 00 if q = 2) with supp b of finite measure for a.e. t G (0,T). Then {9,u) is unique. 

For the proof of Theorem 13.11 only L^-integrability of in time is required. Assuming that 9 
is not only L^-integrable in space uniformly in time, but also L^^'^-integrable uniformly in time we 
obtain the following corollary. 

Corollary 3.3. Let 9(x, t) be a weak solution of the equation liDL^) such that 

9 G L°°((0,r),L2(M'^)) nL«((0,T),LP(M")) withp,qe (2, 00). 

Assume in addition V9 = g + b where g G ^^((0, T), L°^(M")) and b G L^{{0,T), BMO) such that 
supp b has finite measure for a.e. t G (0,T). Then {9,u) is unique. 

The proof of this corollary follows the same outline as the proof of Theorem 1.2 in [16] and we 
refer the interested reader to this paper. 

4. Inviscid limit 

Theorem 4.1. Let 9Q(x,t) be a weak solution of WT^) and let 9f^{x,t) be a family of weak solutions 
of ^DDT^) for some sufficiently small k > 0. Furthermore, assume that 

||^o(-,i)||L2 + \\9^{;t)\\L2 < aoit) and oq G L°°(0,r), 

\\Ooi;t)\\H^ + \\9^{;t)\\H-. <go{t) and go e L\0,T). 

Let s > where q > 2 (s = 00 if q = 2) is such that 9 G L'?((0, T), LP(]R")) for some p G 

(2,;^). IfV9o = g + b with 

\\gi;t)\\L^ + \\bi;t)\\BMO<fo{t) andL%0,T), (25) 
where supp b is of finite measure for a.e. t G (0,T), then 

sup \\{9^,- 9o){-,t)\\L2 ^ 0, asK^O. 

0<t<T 



(24) 
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Proof. Let 0^ and be solutions to problems (|DDEp and (jPEp . respectively. Their difference, 
9 = — Oq, satisfies 

Ot - K{-Aye^ + u^-ve + k - ^^o) • v^o = o, (26) 

where and uq are drift velocities corresponding to 6^ and 6*0- Multiplying by 6 and integrating 
in space using the divergence- free condition on drift velocities, we have 

1 d 



2dt 



dx+ u^- ve^e dx + K K^e^k^e dx = o. 



(27) 



(28) 



Let a = Uf^O and /3 = V^o- We split q = + Or so that |qk| = min(l/K, |q|). Proceeding as in the 
proof of Theorem 13.11 we denote f{t) = 2CCi(||/3i(t)||/,oo + \\(32{t)\\BMo), r{t) = Ci |ar/32| dx 
and x{t) = Ci 11^11^2, and obtain the following differential inequality 

x'{t)<f{t)x{t) (^l + jlnx(t)|+ln(^l + i^^ +r{t) + Kg^{t), 
x(0) = 0. 

To find a good estimate on x{t), we introduce the following equation 

y'{t) = f{t)y{t) (^1 + |lny(t)| +ln (l + ^)) + r{t) + Kgiit), 
y{0) = K 

for some sufficiently small k. From the Osgood existence theorem we know that there exists a 
unique local solution to (129]) . The solution of ([29]) dominates x{t), i.e. 

0<oc{t)<y{t) forte[0,ro], 

where Tq is chosen by similar rules as in the proof of Theorem 13. 1[ From ()29p by Gronwall's 
inequality we have 



(29) 



y(t)<K;exp (^J f{s 



1 + |lny(s)| +ln 1 + 



1 



ds 



1 + llny(T)| +ln 1 + 



dr I ds 



< K exp 1 + 



+ exp 1 + 



In (l + ^)) f{s) ds^ exp f{s)\lny{s)\ ds^ 
In (l + ^)) l^fi^) ds^ exp (^J^ f{s)\ In y{s)\ ds^ ■ {r{s) + gi{s)K) ds. 



(30) 



The condition y{t) > n for sufficiently small n gives | lny(t)\ < — Iuk = ln(l/K). Also let k be small 
enough so that - (l + < thus we can estimate the right hand-side of (|3U|) (up to a constant) 



+ ^ 



{r{s) + go{s)K) ds 



So /(s) ds 



(31) 



K+Cl 



\ar(!^2 \ dxds + k go{s) ds 



Therefore we need an estimate on the part involving a^. As before, we use the fact that with k — ?> 0, 
the support of shrinks. From assumptions on 0^ and we have 



thus by interpolation 



a e L°°((o,r),Li(M")) nLi((o,r),L^(M")), 



a G L5((0,r),LP(M")). 



(32) 
(33) 
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We proceed as in the proof of Theorem 13.11 and, by Holder's inequality, have 

\ar(32\ dx < llar||LP-«ll/?2|lL{p-.)'(3„pp«^) < C||ar H^p-. ||/32 ||iJAfO (34) 

for some sufficiently smah e > (so that p — e > 1). From Chebyschev's inequahty we obtain 

\suppar\ < (K||ar||Lp)^ < («;||a||Lp)^ . (35) 

By Holder's inequality 



l«r-||LP-E = ( / IttrT ) < {supparl^^'''^^ WarWlP (36) 

' supp ar 



Combining (j35|) and (|36|) we obtain 



|ar.||ip-e < np-^\\ar\\^p'' " < KP-E ||a||^p " (37) 



Hence 



I 



ar(32\ dx < KP-'WaWj^p" " \\f32\\BA40 = ""Wall Ip' W Mb pi o- (38) 
Therefore 

Jo Ir '^^^'^ - ''"'''^''((°'^)'-^''(''^"))''^^'''f^'''''~''/''''((0.2'^ ^^^^ 

Thus, if we have (32 £ -^^'^((0, T), BMO) for some s such that g' < s < cx), there exists an e > such 
that above conditions are satisfied and s' = q • We denote 

M = 11/32 ||l=({0,T), BMO)- 

Then, estimates ([30])- ([39]) give 

y{t) < C ds ^ K^C7iM||a||i,((o,T),Lp{M"))) • (40) 

We choose < ti < Tq smah enough so that 1 — 2 J^^ /(s) ds < We get 

y{t) < KkI^- for < t < ti, (41) 

where X is a constant independent of k. We can continue the process on an interval starting at 
t = ti but with the initial condition in ([29]) changed to y{ti) = Kkp-" . Due to integrability of f{t), 
iterating this procedure, we eventually cover the whole interval [0,T]. Letting k — )• we obtain 

sup ||6'k - 6*011^2 ^> 0, as K ^> 0. (42) 

0<t<T 

□ 

Remark 4.2. From the proof of Theorem 14.11 follows that if we assume 

sup \foit) + g^oit)\<N, 

0<t<T 

then we obtain the following explicit rate of convergence 



sup \\e^ - 9o\\l^ < Ck"" 

0<t<T 



The proof follows as in [16 
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Appendix A. Existence of weak solutions to (IDDED 

We sketch the proof of existence of global Leray-Hopf weak solutions of equations (|DDEp . We 
follow the general outline used to construct weak solutions of the Navier-Stokes equations (cf . [20] ) • 
Let (f) G Cq°(IR") be positive, with cp dx = 1. Then cpeix) = e~'^4>{x/e), for e > 0, is a standard 
family of mollifiers. We consider the approximating system 

dtO' + {u' ■ V)0' + K{-Aye' = eA^e' (43) 

with initial data 9q = (p^ * 6q. The drift velocity is obtained from 0^ via a Calderon-Zygmund 
singular integral operator, defined as a Fourier multiplier with a zero-order symbol. Note that 
ll^olli^ — II^o||l2 for any e > 0. 

Let s > § + 1 and fix e > 0. Notice that {-Ay/'^eo G L^{R"-). Since eA^ provides enough 
dissipation, from standard energy estimates it follows that 

sup \\d^{t)\\H^ < C^^n,T,<P,\\eo\\,2' 

te[o,T] 

where Ce^ri,T,0,||6»o||^2 ^ strictly positive constant which is finite for any < T < oo. This a priori 
estimate and a classical Galerkin approach ensure the global existence of a strong H'^ solution to 
(j43p . Moreover, we have the energy estimate 

\\e^{t)\\L2 + K f\\e%s)\\jj,ds<\\eo\\l„ (44) 

Jo 

which is uniform in e > 0. Thus 

r G C7([0,T],L2(IR"))nL2((0,T),F^(IR")). (45) 

This guarantees that, up to a subsequence, 9'' converges weakly-* and weakly to some function 9 
in the spaces L°°((0, T), L2(M")) and L^{{0,T), H"' (W")), respectively. This however is not enough 
to conclude that 9 solves the equation (IDDEp or (IDEp . The classical condition that allows one to 
conclude that 9 is indeed a solution is to establish uniform control over dt9'' in some negative Sobolev 
space. This is accomplished using the energy inequality (I44p . From the Aubin-Lions compactness 
lemma we obtain that 

9' ^9 strongly in L'^{{0,T), Ll^iW^)). 
The strong convergence is necessary in order to pass to the limit in the nonlinear term in the weak 
formulation of (IDDEp . This is done in a standard way using the fact that Calderon-Zygmund 
operators map — )• strongly. 

Acknowledgements. We wish to express gratitude to Susan Friedlander for fruitful discussions and comments on 
early versions of this draft. 
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